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We elaborated the gravitational collapse of a self-gravitating complex charged scalar field in the 
context of the low-energy limit of the string theory, the so-called dilaton gravity. We begin with the 
regular spacetime and follow the evolution through the formation of an apparent horizon and the 
final central singularity. 
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I. INTRODUCTION 

The long-standing prediction of general relativity is the occurrence of spacetime singularity inside black holes. The 
singularity theorems of Penrose and Hawking [1] predict the occurrence of spacetime singularities inside black holes 
under very plausible assumptions, but tell nothing about the geometrical and physical nature and properties of the 
emerging singularities. 

Until recently the only known generic singularity was the Belinsky-Khalatnikov-Lifshitz (BKL) one [2]. According 
to this picture spacetime develops a succession of Kasner epochs in which the axes of contraction and expansion change 
chaotically. This singularity has a strong oscillatory character, which is highly destructive for any physical object. 
Soon after Belinsky et al. [i used a scalar field to have an insight into a cosmological singularity problem of BKL. They 
found that scalar fields destroyed BKL oscillations and singularity became monotonic. Recently, it was shown |4] that 
the general solutions near spacelike singularity in superstring theories and in M-theory (the Einstein-dilaton-p-form 
field) exhibit an oscillatory character of BKL type. 

For a Schwarzschild black hole it has been shown that the asymptotic portion of spacetime near singularity is free 
of aspherical perturbations propagated from the star's surface since the gravitational radiation is infinitely diluted as 
it reaches the singularity. On the other hand, the internal structure of Reissner-Nordstrom (RN) or Kerr black hole 
differs significantly from the above picture. The singularity becomes timelike and both of these spacetimes possess 
Cauchy horizons (null hypersurfaces beyond which predictability breaks down) . In the last few years a new picture of 
the inner structure of black holes was achieved, according to which the Cauchy horizon inside RN or Kerr black hole 
transforms into a null, weak singularity, i.e., an infalling observer hitting this null singularity experiences only a finite 
tidal deformation [5], [5]. The curvature scalars and mass parameter diverge along this Cauchy singularity and this 
phenomenon is known as mass-inflation. The physical mechanism on which the Cauchy horizon singularity is based is 
strictly connected with small perturbations (remnants of gravitational collapse) which are gravitationally blueshifted 
as they propagate in the black hole interior parallel to the Cauchy horizon. For a toy model of spherically symmetric 
charged black hole the main features of singularity at the inner horizon were first deduced analytically from simplified 
examples based on null fiuids f^-[Hj. 

Most of the conclusions supporting the existence of a null weak Cauchy horizon singularity were obtained by means 
of a perturbative analysis. The full non-linear investigations of the inner structure of black holes were given in Ref.[5] 
where the authors revealed that a central spacelike singularity is located deep inside a charged black hole coupled to a 
neutral scalar field. The existence of a null mass- inflation singularity was established by Brady and Smith [HI , in their 
studies of non-linear evolution of a neutral scalar field on a spherical charged black hole. Burko TIT also studied the 
same problem and found a very good agreement between the non-linear numerical analysis and the predictions of the 
perturbative analysis. Expressions for the divergence rate of the blue-shifted factors for that model valid everywhere 
along Cauchy horizon, were given analytically in Ref. 12 . One should have in mind that all these numerical works 
were beginning on RN spacetime and a black hole formation was not calculated. Numerical studies of the spherically 
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symmetric collapse of a massless scalar field in the semiclassical approximation were conducted in Ref. [TH] ■ Piran et al. 
[I4] studied the inner structure of a charged black hole formed during the gravitational collapse of a self-gravitating 
charged scalar field. They started with a regular spacetime and conducted the evolution through the formation of an 
apparent horizon, Cauchy horizon and a final singularity. The results obtained in |14) were confirmed and refined in 
Refs. 

The effect of pair creation in the strong electric fields in a dynamical model of a collapse of the self-gravitating 
electrically charged massless scalar field was elaborated in [T^. The authors studied the discharge below the event 
horizon and its influence on the dynamical formation of the Cauchy horizon. On the other hand, the dynamical 
formation and evaporation of a spherically charged black hole supposed initially to be non-extremal but tending 
towards the extremal black hole and moreover emitting Hawking radiation was investigated in Ref.|18j. Recently, a 
spherically symmetric charged black hole with a complex scalar field, gauge field and renormalized energy momentum 
tensor (in order to take into account the Hawking radiation) was considered in |19| . When the Hawking radiation 
was included it turned out that the inner horizon was separated from the Cauchy one. Studying the neutralization 
of the charged black hole in question it was found that the inner horizon evolved into a spacelike singularity. Using 
the exponentially large number of scalar particles it happened that one could extend investigations inside the inner 
horizon. Recently, the response of the Brans-Dicke field during gravitational collapse of matter was analyzed [20] 
while the internal structure of charged black hole including Hawking radiation and discharge was elaborated in [21] . 
In Ref. [22] the first axisymmetric numerical code testing the gravitational collapse of a complex scalar field was 
presented. Also the non-linear processes with the participation of an exotic scalar field modeled as a free scalar field 
with an opposite sign in the energy-momentum tensor, were considered [23] due to the case when RN black hole was 
irradiated by this kind of matter. 

In this paper we shall consider the implication of superstring gravity for the dynamical collapse of charged complex 
scalar field. The famous Wheeler's dictum that black holes have no hair predicts that the inner structure of the 
black hole will not depend on the collapsing fields. In mathematical formulation this conjecture corresponds to the 
so-called black hole uniqueness theorem [24] (classification of the domains of outer communication of regular black 
hole spacetimes). On the other hand, various aspects of uniqueness theorem for four-dimensional black holes in the 
low-energy string theory was widely treated in Refs. ^25). 

In previous papers |26j we have examined the intermediate and late-time behaviour of matter fields (scalar and 
fcrmions) in the background of dilaton black holes. In some way the present studies will generalize the former ones 
to the more realistic toy model of the dynamical collapse. In what follows we assume that the considered Lagrangian 
for the charged complex scalar field will be coupled to the dilaton via an arbitrary coupling, i.e., e^"'^£(V', -^)^ in 
the string frame. 

The outline of the remainder of the paper is as follows. In Sec. II we derive the equations describing the collapse 
of the charged scalar field in the presence of non-trivial coupling to the dilaton field. Sec. HI was devoted to the 
numerical scheme applied in our investigations. We discussed the numerical algorithm, an adaptive grid used in 
computations and the boundary and initial conditions for the equations of motion for the considered problem. We 
also paid attention to the accuracy of our numerical code. Sec. IV is assigned to the discussion of the obtained results. 
In Sec.V we concluded our researches. 



II. DILATON BLACK HOLE 



In this section our main interests will concentrate on the behaviour of the collapsing complex charged scalar field 
when gravitational interactions take a form typical for the low-energy string theory, the so-called dilaton gravity . To 
take into account the unknown coupling of the dilaton field to the considered charged complex scalar field, we choose 
the action in the form 

j^i?_2(v0)' + e2"^£^ , (1) 
where the Lagrangian C is given by 

C = -\ (V„V + leA^i) (V^,^* - leApr) - F^.F^^r (2) 

The action is written in the string frame but it will be useful to rewrite it in the Einstein frame. In Einstein frame 
the metric is related to the string frame via the conformal transformation of the form provided by the following: 




gap = e "^ffa/?- 



(3) 
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The gravitational part of the action ([T]) appears in the Einstein frame in a more famihar form. Namely, it implies the 
following: 



The equations of motion derived from the variational principle yield 
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' [vpil; + ieAfi^j^ (^V^V* - ieA^i;*^ - ^ae^°"^F^ = 0, 



ieij* ( VtA + ieA^^j ] - ie^ V* - ieA^'il)* 

0, 



= 0, 



(4) 
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where the energy momentum tensor Tpy((/), F, ip, -0*, A) for the fields in the theory under consideration is provided by 
the relation 



In Eq.(lO) by T^itltp^ip* , A) we have denoted the following expression: 

f^,{i;,r,A) = 



ie V '^^V'* + A^ Vf^tP* - ie V* '^^^V' + A^ V^V 



(10) 
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where the explicit form of the Lagrangian Clip , tp* , A) is written as 

C{i;,i;*,A) ^ + ieAp^P^ (^^^/j* - ieA^^^ 

On the other hand, for T^^{F^ (f>) one gets 

1 



T^,u{F, 



2Ff,pF,P ~ -g^,F' ) - 5m- (V0)' + 2V^</)V,<; 



(12) 
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In order to study the gravitational collapse in a spherically symmetric spacetime it will be useful to consider the line 
element written in the double null form |2T 



(14) 



where u, v are advanced and retarded time null coordinates. The null character of the coordinates in question will be 
preserved by the gauge transformation of the form u — f{u) and v — >■ g{v). Using doubly null coordinates enables us 
to begin with the regular initial spacetime at approximately past null infinity, compute the formation of black hole's 
event horizon and then prolong the evolution of the black hole to the central singularity formed during the dynamical 
collapse. 

The assumption of spherical symmetry and the above coordinate choice imply that the only non- vanishing component 
of the C/(l)-gaugc strength field is Fuv or F^^- Consequently, it provides another restriction on the gauge potential. 
Namely, one has to do with A„ or Ay. We can get rid of one of these components of gauge potential by using the 
gauge freedom of the form A^ — ^ + Vud- If one chooses d = J A^dv, then we are left with the only one component 
of the gauge field, which is the function of u and u-coordinates. 

To proceed further, we take into account the w-component of the generalized Einstein-Maxwell equations. It leads 
to the following relation: 

ro„2a4> ^2 A 1 ^2 „24>(a+l) / \ 

(15) 
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Let us define the quantity 



A r-2 



(16) 



just as in Ref415L Q corresponds to the electric charge within the sphere of the radius r{u,v). The above definition 
enables us to separate the second order partial diff^erential equation for into two much simpler first order differential 
equations. We arrive at the following: 



and 



A -^-0 



The equation of motion for dilaton field ^ has the form provided by 



(17) 
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Consequently, the relations for the complex scalar fields are given by 

ic Q (x^ 

r,ui>,v + r yip^u + rijx^uv + ie r Au + ie r y A^ ii H V' = 0, 

r,uVv + r,vi'*u + i~i^*uv - ie r Au ^\ - ie r^^ Au ip* - ^ " ^* = 0. 
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0. (19) 

(20) 
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Combining the adequate components of the Einstein tensor and the stress-energy tensor for the underlying theory we 
obtain the following set of equations: 
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(22) 
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Moreover, we introduce new auxiliary variables written in the form as 

e — ; ^ — ; / — 

a a 

4>, X = (t>^u, 

and the additional quantities provided by the relations as follows: 

M ^ fq + 9P, 
V = gx + fy- 



9 = r^u, 
P = Au, 
y = (l>,v, 



(26) 



(27) 

(28) 
(29) 



Instead of considering two complex fields tp and ip* one can introduce two real fields obeying the relations 4' = 4'i~^i 4^2 
and tp* = — i ip2. On this account, it leads to 



S = SI+IS2, P = Pl+«P2, g = 9l+«g2, 

/i = fj.i + ifi2, = m ^ fq2 + gp2- 



(30) 
(31) 
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Thus, having in mind ah the above, one can rewrite the system of the second order partial differential equations as 
the first order one. By this procedure we get the following system of the first order differential equations: 



PI 
P2 
P3 
P4 
P5 



(He) 



a^u — a c = 0, 
a^y — a d = 0, 
r,u-f = 0, 
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P^iR.) ■ - 91 = 0, 
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P7 : k^u-x = 0, 
P8 : ky-y = 0, 
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(32) 
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Let us introduce some quantities of physical interest. Namely, we define the mass function provided by the relation 



m{u,v) = - 1 + 



4 r.u r. 



1 + 



;2 f9 



(56) 



It represents the Hawking mass, i.e., the mass included in a sphere of the radius r{u,v). Moreover, the Ricci scalar 
has the form as 



R {u, v) 



16x y 
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or it can be rewritten in the form which yields the following: 
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III. NUMERICAL COMPUTATIONS 



A. Numerical algorithm 



The system of equations (32)-(55) in the theory under consideration has to be solved numericahy. In order to find 
the solution one should elaborate an evolution of the quantities d, gi, 52, y^ a, si, S2, k, g, r, Q, j3, /, pi, p2 and 
X. The quantity c does not play the significant role in the process under consideration, so it can be ignored. The 
evolution of the quantities d, qi, q2 and y along u is governed by relations E'i^'^\ "^(fe)' "^(/ni) ^'^'^ D'^'^\ respectively. 
The remaining quantities, a, si, S2, k, g, r, Q, (3, /, pi, p2 and x evolve in turn along w-coordinate according to 
equations P2, P6(fle), P^(im), P^, E2, P4, A/2, Ml, E3'^^\ S^^^^y S"!)^^ and On the other hand, Eqs. PI 

and EA^^') describing the behaviour of c may be discarded and the remaining relations can be used to determine the 
boundary conditions. 

In our studies the numerical algorithm similar to the one proposed in was implemented. The computations 
were carried out on the two-dimensional grid constructed in the (i;u)-plane. In order to obtain a value of a particular 
function at a point (w, it) one should have values of the appropriate functions at points {v — hy, u) and {v, u — hu), 
where /i„ and /i„ are integration steps in v and u directions, respectively. Equations describing the evolution of the 
considered quantities along the coordinates u and v may be symbolically written as: 

fn = F(f,g), g,, = G(f,g). (59) 

In order to get the value of the particular function at a point (v,m), we should find the auxiliary quantities which 
yield 

= f|(.,n-/.„)+^«F(f,g) |(,^_^_^^), (60) 

gg I = g I + Y (g (f, g) I + G (fF, gg) I („ . (61) 

By virtue of the above the final values of the quantities in question are provided by the following: 

f|(.,n) = ^(ff|(„,„)+f|Ku-h„)+'*"F(ff,gg) |(^,^^^, (62) 
g I (v,u) = ^ I (v,u) + g I (v-K^u) + ^-"^ (ff' gg) I ■ (63) 

In the early stages of the calculations the numerical grid is divided evenly, both in v and u directions. On this account, 
at the beginning the quantities and /i„ are equal to each other. 



B. Adaptive mesh refinement 



The coordinates u and v ensure the regular behaviour of all the considered quantities within the domain of integration 
except the vicinity of r = 0. However, during the numerical analysis the considerable difficulties also arise close to 
the event horizon, where function / diverges. A relatively dense numerical grid is necessary in order to satisfactorily 
determine the location of the event horizon and to examine the behaviour of fields inside it, especially for large values 
of the u-coordinate. The efficiency of the calculations suggests using an adaptive grid and performing integration with 
a smaller step in particular regions. On this account the refinement algorithms enable us to make the grid denser 
both in V and u directions and to do the same only along u-coordinate. The first one makes the integration steps /i„ 
and hu smaller on the equal footing as one approaches the event horizon and reaches the large values of v. The other 
one changes only the value of hu- It turned out that the latter refinement of the adaptive grid gives the satisfactory 
results and it is more effective due to the computation time and required computer's memory. Hence, all the results 
presented in our paper will be based on it. 

In order to determine the area of the integration grid, where the grid should be denser, a local error indicator need 
to be used. The aforementioned quantity ought to be bounded with the evolving quantities as well as it should change 
its value significantly in the adequate region. It happened that |15j the function Ar/r along u-coordinate meets our 
requirements. 
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TABLE I: Initial profiles of field functions 



Family 



Profile 



ifD+s) 

ifs) 



p ■ V ■ e 



p- sm^ [-i^fj] ■ cos (Trfyj +2 cos + 5 



C. Boundary and initial conditions 



Having specified the numerical algorithm for the solution of the equations of motion we refine our studies to the 
case of the initial and boundary conditions for the equations in question. The boundary conditions refer to the surface 
u = V, while the initial conditions are formulated along an arbitrarily chosen constant surface u = Ui. A point (0, 0) 
is chosen to be the intersection of these two lines in the (wu)-plane. It also indicates the center of the considered 
spacetime. 

The physical situation we shall take into account will be the gravitational collapse of a spherically symmetric shell 
of infalling complex charged matter. Spacetime of a spherical shell of matter is flat in two regions, i.e., inside the 
shell and at large radii from it. This fact enables us to assume that the line u = v will be not significantly affected 
by the presence of the collapsing shell of matter. On this account the spacetime may be considered as nearly flat and 
electrically neutral in that region. Thus, it gives the following boundary conditions: r — Q = j3 = Q. Consequently, 
equations i^S*^^'^^ reveal that A = along u — v. This fact together with Eqs.PS and P4 provides that / = —g = — f ■ 
Furthermore, the requirement si^r — S2,r = fc,r = = along u = v guarantees that the field functions flatten near 
r = making the numerical analysis possible. Because of the fact that r changes non-linearly along the coordinates 
u and V this boundary condition is implemented using the three-point regressive derivative method with a variable 
step, except the first point, where the Euler's method is used. 

(1 2) 

Combining the relation /^i = /i2 = = along u — v line and Eqs.S'^^'^-j , we obtain that pi — qi, P2 — 92 and 
X — y. On the other hand, the boundary conditions for the quantities evolving along u can be achieved according 
to the algorithm described in the previous section. It can be seen that the auxiliary quantities (61) are taken to be 
equal to the corresponding functions at r = 0. 

The assumption of the fiat geometry in the region, where an observation of the collapse begins justifies the condition 
that d{v,0) = 0. Further, it fixes the remaining freedom in w-coordinate. By virtue of the flatness of the spacetime 
in the vicinity of a surface u = v and the above assumptions we get that a (w, 0) = 1. One should mention that the 
initial conditions ought to include the arbitrary proflles for scalar functions si (v, 0) and S2 (f , 0) describing the real 
and the imaginary parts of the charged scalar field and for the dilaton field k{v,Q). The one-parameter families of the 
initial profiles are listed in Table |Tj where a free family parameter is denoted as p, constants ci and C2 are arbitrarily 
chosen, while Vf equals the final value of w. 6 € [O, ^] is a phase difference determining the amount of the initial 
electric charge ^JJ- The quantities qi (w,0), q2{v,0) and y{v,0) are computed analytically using equations P6(Be}^ 
PQ{im) and P8. The values of the functions g, r, Q, j3, f, pi, p2 and x along the axis u = are obtained using the 
three-point Simpson's method apart from the flrst point, where the Newton's method is implemented. 



D. Numerical tests 



The most straightforward manner of checking the correctness of the numerical code will be a comparison between 
achieved results and an analytical solution of the considered problem. Unfortunately, because of the lack of the 
analytical solution of the problem in question one should apply indirect methods of checking the accuracy of the 
numerical code. 

The first trial will be checking of the convergence of the obtained results. One verifies that the fields under 
consideration converge to some values in the expected region of convergence. From the numerical point of view it 
envisages the fact that the algorithm and its implementation are free of mistakes. To begin with we carried the 
computations not requiring adaptive mesh on four different grids with integration steps equal in both directions. The 
integration step of the particular grid was twice the size of a denser one. The evolving field proflles for arbitrarily 
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chosen u-coordinate are shown in FigjT] We scaled up the vicinities of the cusps, where the differences among profiles 
were most significant. For all the field profiles the very good agreement of an order of O.Of % was achieved. On the 
other hand, the linear convergence of the numerical code is presented in Figj2j The differences between the profiles 
obtained on the two grids with a quotient of the integration steps equal to 2 and their respective doubles are hardly 
distinguishable. The divergence is at most 1%, as may be inferred from the values shown in the magnified areas. 
Moreover, Figj2] gives the clear evidence that the errors become smaller as the grid density increases. 

The next test of our code is to check whether mass (56 1 and charge (161 are conserved in the evolving spacetime. 
The considered fields are scattered by the gravitational and electromagnetic potential barriers as the collapsing 
shell approaches its gravitational radius. Therefore the conservation laws are not satisfied in the entire domain of 
integration. Nevertheless, this effect of the outgoing fluxes of mass and charge is negligible [15] and it has no significant 
influence on the total mass and charge, excluding the area in the vicinity of new forming black hole event horizon. In 
FigjSjthe behaviours of mass and charge for the large value of the advanced time are presented. It turned out that for 
u not exceeding 1, mass and charge are conserved up to within 1.6% and 2.5%, respectively. Further, the inspection 
of FigjS] reveals the deviation of the aforementioned quantities from constancy increases with the advanced time. It 
is caused due to the fact that the reflected waves carry off some mass as well as charge. Because of the fact that the 
total mass contains also the energy momentum of gravitational field which is not taken into account in the Noether 
current bounded with the energy momentum tensor, the subject of the mass conservation is more subtle [15j . 

The last test of the accuracy of our code consists of the analysis of the simplified versions of the problem in 
question. Namely, in Fig|4]we depicted the outgoing null rays in (ru)-plane for the spacetime containing black hole 
stemmed from the gravitational collapse of the neutral scalar field. The situation corresponds to setting a = e = 
in the equations of motion and eliminating electrically charged scalar field by putting ■0 = 0. The value of the free 
parameter is equal to 0.125. On the other hand, Figj5] illustrates the outgoing null rays for the spacetime of a black 
hole emerging due to the gravitational collapse of an electrically charged scalar fields. Here, we put a = and get rid 
of the dilaton field by setting = 0. The value of free parameter was taken to be 0.5. The inspection of Figj5] shows 
the formation of the black hole event horizon and a Cauchy horizon at asymptotically large v. The structures of the 
spacetimes emerging during collapses presented in Figs|4]and[5]are in a perfect agreement with the results published 
in Refs.tlIl[Tl[Tn]. 



IV. RESULTS 



In our numerical studies we have used the one-parameter families of initial profiles referring to the real and imaginary 
parts of the electrically charged scalar field and dilaton field. The results do not depend on the type of the family 
of the initial profiles as well as on family constants. Hence their choice is unrestricted. Moreover, no significant 
dependence on the electric coupling constant was observed. On this account the electric coupling constant was put 
e = 0.5, in all the calculations. All the results in the present section were obtained using profiles [fo+s) with values 
of the family constants equal respectively to: ci^si = 0.75, C2,si = 0.15, ci^sj = 1.55, C2,S2 — 0.17 and ci_fe = 1.3, 
C2,fc = 0.21. Subscripts Si, S2 and k refer to the real and imaginary parts of the electrically charged fields and dilaton 
field, respectively. The family parameters , and pk are taken to be equal and they are denoted by p. 

The critical phenomena in gravitational physics turned out to be one of the key thought experiments in the studies 
of black hole formation (see Ref.[28ji29j and references therein). There is a critical value of the parameter, denoted 
by p*, below which the spacetime is non-singular and does not contain a black hole. For values exceeding p* there is 
a black hole in the spacetime, which means that it is singular. These phenomena are referred to as subcritical and 
supercritical, respectively. It was revealed in Ref . ^30j that for any fixed radius observer, as we take the limit to the 
critical value of the parameter, the sphere of the influence of the diminishing black hole mass shrinked to zero. On the 
other hand, the resulting spacetime converges pointwise to Minkowski spacetime at r > 0. Moreover, the convergence 
is not even uniform. 

In what follows we shall study the gravitational collapse of a self-interacting complex charged scalar fleld in the 
dilaton gravity. We take into account models with different values of coupling constant a. To begin with one considers 
first a subcritical evolution. In Figj6] and [s] we depicted the radial function r{u,v) as a function of the ingoing null 
coordinate v along a sequence of the outgoing null rays (u = const.) All the outgoing null rays originate from the 
nonsingular axis r — 0. Just we begin our evolution with a regular spacetime along u — 0. In Figs|6] and |8] we 
plotted two cases, first one for a slightly curved spacetime when p ^ p* and the other for almost critical one, for 
which p ^ p* . We take into account two values of free parameter p = 0.05 and p = 0.0517693, for the case when the 
coupling constant a = 1 and e = 0.5, and p = 0.0485 and p = 0.0503946, when the value of coupling constant a = — 1 
and e = 0.5. In both figures we observe a flat region corresponding to small values of advanced and retarded times 
in which one has that r{u,v) is proportional to u and v. With the passage of time the curvature of the spacetime 
becomes considerable and we have no longer proportionality between r and {u,v). As was expected the spacetime 
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under consideration curved earlier and more significantly for larger value of the parameter p. In Fig|6]we observe the 
intersections oi u = const, lines for larger values of r. 

Figsj?] and [9] display Kruskal diagrams {r{u,v),v) for the spacetimes under consideration. As in the previous case 
we start the evolution from a regular spacetime along u — 0. We studied two cases of singular spacetime, one with 
small black hole when p* and the other one, with big black hole for which p^ p* . The values of the p parameter 
are 0.0525 and 0.1, respectively. FigjT] describes the birth of black holes for the coupling constant a — 1 and e — 0.5, 
while Figj9]is performed for the case a = — 1 and e = 0.5. In both figures one can distinguish between two types of 
outgoing null rays in the (ru)-plane. Namely, the outermost null rays escaping to the future null infinity (for small 
value of retarded time) and the innermost null rays which stem from the nonsingular axis r = and terminate at the 
singular part of the hypersurface r = 0. Their evolution is described by finite value of v. Contrary to Refs.[14l [15] 
we did not find in our numerical simulations intermediate outgoing null rays approaching a fixed radius at late times 
(when V — >■ 00). Just, the evolution in question resembles formation of Schwarzschild black hole, rather than RN one 
with Cauchy horizon. 

In order to better understand the causal structure of the dynamical spacetimes in question we shall proceed to 
perform the Penrose diagrams (the dependence {u,v) along r = const.). In Figs 10 and 12 we present the results 
for the slightly curved and almost critical non-singular spacetimes for Q; = lora = — 1 and e = 0.5. The other 
parameters are as in Fig|6] and [8] respectively. For both of them the outermost contour line corresponding to r = is 



a non-singular straight line. On the other hand, in Figs 11 and 13 we presented lines of constant r in the (t;u)-plane 



for small and large black hole emerging from the gravitational collapse. Fig |ll| was performed for a = 1, while in 
Fig 13 we put the coupling constant a — —1. In both cases one has e — 0.5. The outermost thick line is equivalent to 
r = 0. Contrary to the previous cases, in the spacetime of dynamical formation of black hole one has a straight line 
(u ~ v) in the left section. This behaviour corresponds to the non-singular axis. On the other hand, the right part 
corresponds to the central singularity at r = 0. Because of the fact that r^v < along the latter section, we have to 
do with the spacelike singularity. Two types of horizons are present in the singular spacetimes. The apparent horizon 
is represented by the contour r ,j ~ 0, while the event horizon is provided by the line of constant u and characterized 
by f y = when w — >■ 00. The dynamical character of the emerging spacetimes is refiected in relative positions of the 
aforementioned horizons. They do not coincide in the early stages of the evolution and in the end, they approach the 
same value of u = const, as v tends to infinity. 

The next object of our interest was an influence of dilatonic coupling constant on the evolution described by the 
considered equations of motion. The set of Penrose diagrams (representing lines r = const, in (uu)-plane) for different 
values of a is shown in Fig |14[ The values of dilatonic coupling constants were arbitrarily chosen to be ±0.5, ±1 
and ±1.5. The value of the free family parameter was taken a.s p = 0.075. We conclude that the value of a exerts 
no qualitative influence on the structure of spacetime. The most striking effect is connected with the moment of the 
horizon's formation. In general, for bigger absolute values of a the horizon forms at earlier advanced times. Although 
the effect is far more noticeable for positive values of dilatonic coupling constant, it is also present for a not exceeding 
zero. 

We also examined mass of a black hole emerging from the gravitational collapse in question as a function of v- 
coordinate along the apparent horizon for different values of coupling constant a. The relations are depicted in Figs|15| 



and 16 Mass of a black hole, denoted by A/, is Hawking mass (56) calculated along apparent horizon. It turned out 



that for large values of retarded time black hole mass tends to a constant value. For positive values of a we observe 
that the bigger coupling constant is the bigger Hawking mass one gets. Moreover, the dependence of M on a-coupling 
constant is linear. In the case of a less than zero, its influence on the asymptotic value of Hawking mass is not so 
st raight forward . 

Now we proceed to study different features of the collapsing spacetimes as functions of null coordinates. Namely, 
we shall concentrate on Hawking mass (56), Ricci scalar (57), metric coefficient g^v and w-derivative of r. In our 
considerations we put p = 0.1, e = 0.5 and considered two distinct types of evolutions for dilatonic coupling constants 
equal to a = ±1. We took into account foliations of spacetime in both v- and u-directions and studied ingoing null 
rays terminating at a non-singular part of r = 0, outside event horizon and within it, as well as at singular r = 0. We 
paid attention to outgoing null rays outside, inside and exactly along the event horizon of the emerging black hole. 
For both values of a we examine ingoing null ray v — 0.75 lying entirely outside the event horizon, hence terminating 
at the non-singular part of r = 0. Next, we consider a null ray situated along v = 1 which crosses the event horizon, 
but also ends at regular r = 0. The other ingoing null rays are v — 1.25 for a — —1 and v — 1.2 for a — 1. They have 
the same characteristic as the previous one, but they are situated much closer to the cusp on Penrose diagram of the 
spacetime. We also consider two null rays terminating at central singularity which lie along v = 1.75 and v = 2. As 
far as the outgoing null rays are concerned, one examined the outgoing null ray u = 0.5 situated outside event horizon 
for a = —1 and a = 1. The event horizon is situated along u = 0.8557 for a = — 1 and along u = 0.6663 for a = 1, 
so these null rays were the next objects of our studies. We also paid attention to u = const, lying beyond the event 
horizon. They are situated along u — 0.86, u — 0.95 for a — —1 and along u = 0.67, u — 0.75 for a = 1. 
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Hawking mass as a function of u along null rays in question for a = —1 and a = 1 is shown in Fig 17 For both 



values of a, in the case of ingoing null rays terminating at non-singular part of r = 0, Hawking mass tends to zero 
as the retarded time increases. Its behaviour along v — const, is evidently different for ingoing null rays hitting the 
central singularity. For a — — I Hawking mass is constant until it reaches the vicinity of r = 0, where it grows rapidly. 
The situation is quite different for the coupling constant a = 1. In this case Hawking mass grows continuously along 
V — const, and its increase is not so sudden close to r = 0. 

Hawking mass as a function of v along the aforementioned null rays for a = —1 and a = 1 is shown in Fig |18| 
It grows during the dynamic evolution and reaches a constant value when the spacetime becomes static, which is 
indicated by the coalescence of apparent and event horizons. In the early stages of the collapse, for both a, the values 
of Hawking mass are bigger for small retarded times. For a = —1 at some z;- interval this relationship disappears and 
Hawking mass randomly reaches actually the same constant value for all outgoing null rays. The interval is shown in 
the magnified area in Fig |18^ . On the contrary, for a — 1 there is one point, at which lines indicating Hawking mass 
as a function of v cross. The relationship mentioned above changes its character and values of Hawking mass are now 
bigger for large retarded times. The point and its vicinity are depicted in the magnified area in FigfT8|3. Moreover, 
for a = 1, the constant values of Hawking mass for large advanced times are different for each u — const. 



Ricci scalar as a function of u and v along the considered null rays is shown in Figs 19 and 20 In flat regions of 



the spacetime its values are practically equal to zero. For a — —\ and a ~ 1 Ricci scalar is finite along ingoing null 
rays terminating at the non-singular r = and along outgoing null rays escaping to infinity. It is divergent along null 
rays ending at the singular part of r = 0. 



On the other hand, the metric coefficient Quv as a function of u and v along null rays in question is shown in Figs 21 



and 22 For both values of a it is constant in flat regions of spacetime. In curved areas it is slightly increasing with 
the retarded time for ingoing null rays terminating at a regular part of r = and considerably decreasing in the other 
case, guv has a peak near the value of the advanced time, where the apparent horizon appears and the spacetime 
becomes singular. The peaks for both values of a are shown in magnified areas of the respective figures. Then, 
decreases with the advanced time along all outgoing null rays. The changes are rather small outside the apparent 
horizon and become large beyond it. 

The derivative of r with respect to w as a function of u along null rays in question is shown in Fig |23| It is constant 
or slightly decreasing for ingoing null rays lying outside the apparent horizon and displays a strong variability along 
V = const., which crosses the apparent horizon. On the other hand, the derivative of r with respect to f as a function 
of V along the examined null rays is shown in Fig |24[ It is constant in the w-interval corresponding to the non-singular 
part of r = 0. In the region, where the line indicating a central singularity on Penrose diagram is the steepest 
decreases rapidly, but not uniformly - it in turn displays quick and slow changes. In the area, where the singular r = 
lies almost along a constant u the derivative of r with respect to v changes only a little. It increases for outgoing null 
rays outside the event horizon, is constant along it and decreases within it. 



V. CONCLUSIONS 



In our paper we considered a dynamical collapse of a charged complex scalar field in the low-energy limit of the 
string theory, the so-called dilaton gravity. In order to solve the complicated equations of motion one used the double- 
null coordinates which enabled us to begin with the regular initial spacetime at null infinity, compute the formation 
of a black hole and broaden our inspection to the singularity formed during the dynamical collapse in question. We 
have formulated the problem under consideration as a system of the first order partial differential equations with the 
adequate boundary and regularity conditions. 

We begin our numerical studies with a subcritical evolution. One observes that there is a flat region in (rv)- 
plane corresponding to small values of advanced and retarded coordinates where r(u,v) is proportional to u and 
w-coordinates. Then, the curvature of the spacetime comes into being, destroying this behaviour. The spacetime 
curved earlier and more significantly for the larger values of p parameter. Next we proceed to analyze the Kruskal 
diagrams (r, w) along u = const, for the considered spacetimes. We found that the outermost null rays escaped to 
the future null infinity, while the innermost originating from the non-singular axis r — Q terminated at the singular 
part of the hypersurface in question. We did not observe formation of the inner horizon. The causal structure of 
the dynamical collapse was elaborated by means of Penrose diagrams, namely we analyzed (w, M)-dependence along 
r — const. We noticed that the left part of the diagram corresponded to the non-singular axis r — Q and the right 
section of it was bounded with the central singularity. We have indicated two types of the horizons, i.e., the apparent 
one represented by the contour = and the event horizon forming when v ^ oo. They did not coincide at the 
early stages of the dynamical collapse but in the end they approached the same value of M-coordinate as, v co. 
Penrose diagrams for various values of a coupling constant were also found. 

On the other hand, we considered Hawking mass as a function of ?;-coordinate along different null rays. It turned out 



11 



that it grew during the dynamic collapse and reached the constant value for a static spacctime, where the coalescence 
of the apparent and event horizon took place. We also analyzed the behaviour of Hawking mass as a function of 
ingoing null rays terminating at a non-singular part of r = 0. It was revealed that it tended to zero with the passage 
of the retarded time. Hawking mass behaviour, along v = const, for ingoing null rays hitting the central singularity, 
was different for the different values of coupling constant a. 

As far as the Ricci scalar and metric tensor coefficient Quv are concerned their behaviour is more or less similar for 
the positive and negative a coupling constant. 

To conclude one remarks that the dynamical collapse of a complex charged scalar field in the realm of dilaton 
gravity resembles the collapse leading to the Schwarzschild black hole rather than the collapse of charged field in 
Einstein-Maxwell theory. Though, when we check our code and put dilaton field and coupling constant equal to zero 
we get the behaviour leading to black hole with an a Cauchy horizon, during dynamical collapse of charged scalar 
field in dilaton theory we do not find the inner black hole horizon. 

Acknowledgments 

AB was supported by Human Capital Programme of European Social Fund sponsored by European Union. 



[1] S.W. Hawking and G. F.R.Ellis, The Large Scale Structure of Space-time (Cambridge University Press, Cambridge, England, 
1973). 

[2] V.A.Belinsky, I.M.Khalatnikov, and E.M.Lifshitz, Adv. Phys. 19, 525 (1970). 
[3] V.A.Belinsky and I.M.Khalatnikov, Zh. Eksp. Toor. Fiz. 74, 3 (1978). 
[4] T.Damour and M.Hcnncaux, Phys. Lett. B 488, 108 (2000), 

T.Damour and M.Hcnncaux, Phys. Rev. Lett. 85, 920 (2000). 
[5] A.Ori, Phys. Rev. Lett. 67, 789 (1991), 

A.Ori, Phys. Rev. Lett. 68, 2117 (1992), 

L.M.Burko and A.Ori, Phys. Rev. Lett. 74, 1064 (1995), 

A.Ori and D.Gorbonos, J. Math. Phys. 48, 092502 (2007), 

D. Gorbonos and G.Wolansky J. Math. Phys. 48, 092503 (2007). 
[6] P.R.Brady and J.D.Smith, Phys. Rev. Lett. 75, 1256 (1995). 

[7] W.A.Hiscock, Phys. Lett. A 83, 110 (1981). 

[8] E.Poisson and W.lsrael, Phys. Rev. D 41, 1796 (1990), 

E. Poisson and W.lsrael, Phys. Rev. Lett. 63, 1663 (1989). 

[9] M.L.Gnedin and N.Y.Gnedin, Class. Quantum Grav. 10, 1083 (1993). 
[10] L.M.Burko, Phys. Rev. Lett. 79, 4958 (1997). 

[11] R.S.Hamadc and ,J.M. Stewart, Class. Quantum Grav. 13, 497 (1996). 
[12] L.M.Burko and A.Ori, Phys. Rev. D 57, R7084 (1998). 
[13] S.Ayal and T.Piran, Phys. Rev. D 56, 4768 (1997). 
[14] S.Hod and T.Piran, Phys. Rev. Lett. 81, 1554 (1998), 
S.Hod and T.Piran, Gen. Rel. Grav. 30, 1555 (1998). 
[15] Y.Oren and T.Piran, Phys. Rev. D 68, 044013 (2003). 
[16] J.Hansen, A.Khokhlov, and I.Novikov, Phys. Rev. D 71, 064013 (2005). 
[17] E.Sorkin and T.Piran, Phys. Rev. D 63, 084006 (2001). 
[18] E.Sorkin and T.Piran, Phys. Rev. D 63, 124024 (2001). 

[19] S.E.Hong, D.Hwang, E.D.Stewart, and D.Yeom, Class. Quantum Grav. 27, 045014 (2010). 
[20] D.Hwang and D.Yeom, Class. Quantum Grav. 27, 205002 (2010). 

[21] D.Hwang and D.Yeom, Internal Structure of Charged Black Holes, gr-qc 1010.2585 (2010). 
[22] E.Sorkin, Phys. Rev. D 81, 084062 (2010). 

[23] A.Doroshkevich, J.Hansen, D.Novikov, I.Novikov, D.Park, and A.Shatskiy, Phys. Rev. D 81, 124011 (2010). 
[24] M.Heusler, Black Holes Uniqueness Theorems (Cambridge University Press, Cambridge, England, 1996), 

K.S.Thorne, Black Holes and Time Warps (W.W.Norton and Company, New York, 1994). 
[25] A.K.M.Masood-ul-Alam, Class. Quantum Grav. 14, 2649 (1993), 

M.Rogatko, Class. Quantum Grav. 14, 2425 (1997), 

M.Rogatko, ib'id. 19, 875 (2002), 

M.Mars and W.Simon, Adv. Theor. Math. Phys. 6, 279 (2003), 
M.Rogatko, Phys. Rev. D 58, 044011 (1998), 
M.Rogatko, ibid. 59, 104010 (1999), 

M.Rogatko, ibid. 82, 044017 (2010), 
S.Yazadijev, tbid. 82, 124050 (2010). 



12 



[26] R.Modcrski and M.Rogatko, Phys. Rev. D 63, 084014 (2001), 
R.Modcrski and M.Rogatko, 'ibid. 64, 044024 (2001), 
R.Modcrski and M.Rogatko, ibid. 72, 044027 (2005), 
M.Rogatko, ibid. 75, 104006 (2007), 
R.Moderski and M.Rogatko, ibid. 77, 124007 (2008), 
G.Gibbons and M.Rogatko, ibid. 77, 041034 (2008), 
M.Gozdz, L.Nakonieczny, and M.Rogatko, tbid. 81, 1040027 (2010). 

[27] D.Christodoulou, Comm. Pure Appl. Math. 46, 1131 (1993). 

[28] M.W.Choptuik, Phys. Rev. Lett. 70, 9 (1993). 

[29] C.Gundlach, Phys. Reports 376, 339 (2003). 

[30] A.V.Prolov and U.Pen, Phys. Rev. D 68, 124024 (2003). 



13 




FIG. 1: Plots of the field profiles: (a) real and (b) imaginary parts of electrically charged scalar field, (c) dilaton field, along 
u = 0.7 taken for the different integration steps (multiples of /i = 10~*). The evolution was monitored for a Gaussian initial 
pulse {fn+s) with the parameter p = 0.005 and constants (ci,C2) equal to (0.75,0.15), (1.55,0.17), (1.3,0.21), respectively. 
The values of coupling constants are a = 1 and e = 0.5. 
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FIG. 2: The differences between the field profiles: (a) real and (b) imaginary parts of electrically charged scalar field, (c) dilaton 
field and their corresponding doubles for evolutions elaborated in Figjl] The subscripts of field functions denote integration 
steps (multiples h = 10~*). 
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FIG. 3: Mass and charge as functions of u-coordinatc along v = 5 outside event horizon for: (a) a = —1 and (b) a = 1. Initial 
profile of the dilaton field was a Gaussian pulse (fo+s) with constants (ci,C2) equal to (0.17,0.21). The charged scalar field 
was initially represented by the profile (/s) with constants Vf = 5 and <5 = f ■ The value of electric coupling constant is e = 0.5, 
the free parameter is taken to be p = 0.075. 
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FIG. 4: Outgoing null rays in (ru)-planc for the spacctimc containing a black hole left after gravitational collapse of a neutral 
scalar field. The initial profile belongs to the family (fo+s) with constants (01,02) equal to (1.3,0.21). The value of the free 
parameter is p = 0.125. 
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FIG. 5: Outgoing null rays in (r'(j)-planc for the spacctirnc containing a black hole left after gravitational collapse of an 
electrically charged scalar field. The initial profile belongs to the family (fs) with constants equal to Vf = 7.5 and 5 = f . The 
value of electric coupling constant is e = 0.5 and the value of free parameter is p = 0.5. 
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FIG. 6: Outgoing null rays in (r7j)-planc for: (a) slightly curved {p <C p*) and (b) almost critical {p ^ p*) non-singular 
spacetime left after the evolution of the discussed fields. The values of the free parameters are p = 0.05 and p = 0.0517693, 
respectively. The coupling constants are a = 1 and e = 0.5. 
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FIG. 7: Outgoing null rays in (r«)-plane for spacetime containing: (a) a small black hole (p ^ p*) and (b) a large black 
hole (p ^ p*) emerged from the gravitational collapse under consideration. Formation of the apparent and event horizons is 
depicted. The values of the free parameters are p = 0.0525 and p — 0.1, respectively. The coupling constants a and e are the 
same as in Fig|6] 
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FIG. 8: Outgoing null rays in (rw)-plane for: (a) slightly curved and (b) almost critical non-singular spacetime left after the 
evolution of the fields in question. The free paranaeters are taken to be p = 0.0485 and p = 0.0503946, respectively. The 
coupling constants are a = —1 and e = 0.5. 
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FIG. 9: Outgoing null rays in (rw)-plane for spacetime containing: (a) a small and (b) a large black hole emerged from the 
gravitational collapse under consideration. Formation of an apparent and an event horizons is depicted. The values of the free 
parameters are: p — 0.0525 and p = 0.1. The values of the coupling constants are the same as in Fig[8] 
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FIG. 11: Lines of constant r in (iiu)-plane for spacetime containing: (a) a small and (b) a large black hole emerged from the 
considered gravitational collapse. Formation of an apparent and an event horizon and the position of singularity are shown. 
The values of the free parameters and the coupling constants are the same as in Fig[7| 
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FIG. 13: Lines of constant r in (iiu)-plane for spacetime containing: (a) a small and (b) a large black hole emerged from the 
gravitational collapse under consideration. Formation of an apparent and an event horizons as well as the position of singularity 
are shown. The values of the free parameters and the coupling constants are the same as in Fig|9] 
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FIG. 16: The influence of positive values of coupling constant a on mass of emerging black hole: (a) mass as a function of 
v-coordinate along the apparent horizon and (b) an asymptotic value of black hole mass (for u = 3) as a function of coupling 
constant a. 
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FIG. 17: Hawking mass as a function of M-coordinate along ingoing null rays for: (a) a = —1 and (b) a = 1. 
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FIG. 18: Hawking mass as a function of u-coordinate along outgoing null rays for: (a) a = — 1 and (b) a = 1. 
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FIG. 19: Ricci scalar as a function of M-coordinate along ingoing null rays for: (a) a = —1 and (b) a = 1. 



32 



SQ 



-90 



-10O 



-150 



-200 



-2S0 




T" 





w ■ 




m'W^ — 

i 


■ fi 





- i 



' I 




— 



-so 



-100 



-ISO 



-200 



-250 




U-O.S • 

u-o.sseG 

U-0.67 — 





I _ 



If u 

if V 



4 



2j5 



□.5 



1 .5 



2.a 



FIG. 20: Ricci scalar as a function of v-coordinate along outgoing null rays for: (a) a = — 1 and (b) a = 1. 
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FIG. 21: Metric coefficient as a function of u-coordinate along ingoing null rays for: (a) a = —1 and (b) a = 1. 



36 



■1v 



-14 




I 



T 1 r 




J I I I L 



U-Cl.S 
U-0.3557 
U-Q.SS 
U-D.95 



□ .5 1 1.5 



2 2.a 3 3.5 4 4.5 5 
V 



0.6 



0.4 



0.2 



J"- -O.E 



-0.4 



-0.6 



-o.a 



-I r 



T 1 1 1 r 




FIG. 24: Quantity r,^ as a function of v-coordinate along outgoing null rays for: (a) a = — 1 and (b) a = 1. 



